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t^. ' 1 Introduction 

O ■ 

The Classical Analytical Mechanics is an old and well-established part of both mathematics and physics. 
Nevertheless many people still look for the best mathematical tools in describing various aspects 
of mechanical systems. A use of Lie algebroids and Lie groupoids for describing some systems of 
the Classical Mechanics was proposed by P. Libermann [S3] and A. Weinstein [?5] more than ten 
years ago. This turned out to be a very fruitful idea and since then much work has been done (e.g. 
[H Q~2] HH HH HH1 132) making use of Lie algebroids in various aspects of Classical Mechanics and 
Classical Field Theory. The need of extending the geometrical tools of the Lagrangian formalism from 
just tangent bundles to Lie algebroids is justified by the fact that reductions usually move us out of 
the environment of the tangent bundles (think on the rigid body). It is similar to the better-known 
situation of passing from the symplectic to the Poisson structures in the Hamiltonian formalism. 

In the paper [T2] it was observed that, following some ideas of W. M. Tulczyjew and using general 
algebroids instead of just Lie algebroids, one can describe a larger class of systems in a simple and 
elegant way, both in the Lagrangian and in the Hamiltonian formulation. Moreover, the proposed 
geometric picture does not require considering prolongations of Lie algebroids we start with, as it was 
in the case of previous approaches known in the literature. A further paper 13J, was devoted, in turn, 
to the construction of Euler-Lagrange equations in the affine setting of so called special affgebroids 
which is particularly suitable for time-dependent systems. 

In the present paper we concentrate on variational calculus and constraints in the algebroid setting. 
We work with a general algebroid, defined in [17] as a double vector bundle morphism 

(1.1) e:T*E~^TE* 
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covering the identity on E* . Here r : E — > M is a vector bundle playing the role of kinematic 
configurations. To some extent then, our paper can be understood as a natural generalization of [33] , 
where a variational calculus on Lie algebroids has been developed according to the original ideas of 
A. Weinstein [45], and of [5] [18], where constraints on Lie algebroids have been considered. On the 
other hand, our approach is definitely different from the approaches known in the literature, even when 
the equations we obtain cover the corresponding Euler-Lagrange equations in the Lie algebroid case. 
This is mainly because we adapt the framework of the Tulczyjew triple [381 1391 141] , working simply 
with the morphism (jl.ip rather than following the Klein's method [21] generalized to Lie algebroids, 
in which the bundles tangent to E and E* are replaced by the prolongations of E with respect to 
the vector bundle projections r: E — > M and ir: E* — » M. This, in our opinion, simplifies the whole 
formalism substantially. 

To define a variational problem on an algebroid we have to specify a manifold M. of paths whose 
tangent space T M. represents all possible variations and an action functional W on M.. Then we 
have to choose a submanifold M of admissible paths and a set (generalized distribution) V C TM. 
of admissible variations of admissible paths. In [33j admissible variations are constructed out of 
homotopies of admissible paths as defined in [8]. For general algebroids we need different way of 
constructing admissible variations, since we have to accept the fact that they are not tangent to 
the submanifold of admissible paths in general. Therefore we construct admissible variations for an 
admissible path 7 in E out of vertical variations of 7 in E, i.e. out of vertical vector fields along 7. 
Note that the variations are defined in E (which is TM in the standard variational calculus), not in 
M. This is because the variational calculus on algebroids leads to first-order differential equations in 
E rather than to second-order equations in M. This is only the case of the canonical Lie algebroid 
E = TM when paths in M are in one to one correspondence with admissible paths in E, this time 
- just tangent prolongations of paths in M , and admissible variations are tangent prolongations of 
variations of paths in M. For a general algebroid the admissible variations are constructed from the 
vertical ones by means of the double vector bundle relation k = n e : TE — \>TE which is dual to 
the morphism e. Of course, for Lie algebroids our admissible paths coincide with the infinitesimal 
homotopies of admissible paths associated with the lifts of time-dependent sections, as they appear in 
33 , 8 J . We prefer a more fundamental approach which uses n e to produce admissible variations out of 
the vertical ones, instead of lifting whole sections extending paths in E and showing that the result does 
not depend on the extension. In the case of E = TM the mapping e defining an algebroid structure 
is the inverse to the Tulczyjew isomorphism olm '■ TT*M — > T*TM. The relation k £ is in this case the 
well-known canonical flip km ■ TTM — > TTM. Our construction is especially convenient in the case 
of non-holonomic constraints where variations are not tangent to the submanifold of constraints. 

It is clear from our variational picture that putting constraints must result in defining a subset of T>. 
In the case of a general algebroid E our classification of the constraints is based on the way in which the 
constrained admissible variations are constructed. According to the tradition we call them: vakonomic, 
non-holonomic, and holonomic constraints. Starting from a subset S of E, classically understood 
as a geometric constraint for velocities, we have at least two natural possibilities of constructing a 
constraint in admissible variations: one is to consider only admissible variations which are tangent to 
S (vakonomic constraint), the other - to consider only admissible variations coming from those vertical 
ones which are tangent to S (nonholonomic constraint). Note that our approach allows to understand 
nonholonomic constraint as a constrained variational problem, contrary to the commonly accepted 
conviction. A nonholonomic constraint is called holonomic if the constrained admissible variations are 
tangent to S (are vakonomic). Sometimes it is hard to decide without making an experiment which 
method should be used to describe the real behavior of the system. 

For all types of constraints we construct analogs of the Euler-Lagrange equation for systems that 
are subject to those three types of constraints in variational way. Note however that the corresponding 
equations describe "regular" solutions rather than a general solution of the variational problem. Ad- 
ditionally, like for non-constrained cases in |12j , we derive the equations purely geometrically, without 
referring to the variational calculus. 

The literature concerning constraints in Variational Calculus is so extensive that there it is impossi- 
ble to cite it in a complete way. We decided to list among references only papers dealing actually with 
Lie algebroids or being direct inspiration for the framework we propose. Let us also make it clear that 
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we see the meaning of the present paper not only as a generalization of formalisms of Classical Me- 
chanics. Working with the case of a general algebroid forced us to propose a geometric approach which 
seems to be new and illustrative even when applied to very classical situations. The main observation 
is that an algebroid structure is a crucial geometric ingredient in constructing the dynamics of the 
system. It tells us not only the configurations, velocities and inner degrees of freedom, but it contains 
the information on how the admissible variations should be produced from a simple geometric model 
of variations of paths in a vector bundle - the vertical ones. This structure is encoded in a single map 
respecting double vector bundle structures. The brackets and the Jacobi identity are therefore 
proven to play a minor role. The Jacobi identity for an algebroid bracket ensures some integrability 
conditions that allow us to integrate the Lie algebroid into an (at least local) Lie groupoid (see [8]), but 
which is irrelevant for the possibility of constructing Euler-Lagrange equations. Fixing this geometric 
setting for our system, it is then the Lagrangian function which produces a concrete dynamics out of 
these data. However, we would like to stress that regularity of the Lagrangian is completely irrelevant 
for our picture. The general method of constructing dynamics out of the Lagrangian works for all 
Lagrangians, singular or not. The difficulty with singular Lagrangians is that the dynamics we obtain 
is really implicit and complicated. In other words, difficulty with singular Lagrangians lies in difficulty 
in solving equations, not in the geometric construction of the equations themselves. 

Finally, if the variational calculus is concerned, only admissible paths come to the play. This 
is because we work on the bundle E of kinematic configurations and considering only admissible 
paths corresponds, classically, to work with paths in the manifold M of position configurations lifted 
canonically to the paths in TM. The geometrical model of (infinitesimal) variations of an admissible 
path 7 is to consider vertical vector fields along 7 : [£ , h] — ► E. Now, the true (mechanical) admissible 
variations are vector fields along 7 constructed from the vertical ones out of the algebroid structure k. 
This is how the algebroid structure comes to the variational picture. Note that the role of the (Lie) 
algebroid structure in the classical setting is usually overlooked, since it is hidden behind structures of 
the tangent and cotangent bundles which are viewed as a natural part of the theory. 

The paper is organized as follows. In Section 2 we set up the notation and we recall the notion of 
general algebroid as a double vector bundle morphism. Then we introduce the relation k that is used 
for defining admissible variations. In Section 3 we discuss the Lagrange formalism without constraints 
on general algebroid. Then we pass in Section 4 to the variational calculus. We derive the variation of 
the Lagrangian and Euler-Lagrange equations. The final section is devoted to constraints. Geometric 
constraints as subsets S C E give rise to variational constraints which are classified in pure geometrical 
terms as vakonomic, nonholonomic, or holonomic. We derive constrained equations using variational 
motivations and give them pure geometric interpretations. 

2 Lie algebroids as double vector bundle morphisms 

We start with introducing some notation. 

Let M be a smooth manifold and let (x a ), a = 1, . . . , n, be a coordinate system in M. We denote 
by tm ■ TM — > M the tangent vector bundle and by ttm ■ T*M — > M the cotangent vector bundle. We 
have the induced (adapted) coordinate systems (x a ,x h ) in TM and (x a ,pb) in T*M, Let r: E — > M 
be a vector bundle and let tt: E* — > M be the dual bundle. Let (ei, . . . , e m ) be a basis of local sections 
of t; E — > M and let (el, ... , e™) be the dual basis of local sections of tt: E* — > M. We have the 
induced coordinate systems: 



where the linear functions t(e) are given by the canonical pairing t(e)(v x ) — (e(x),v x ). Thus we have 
local coordinates 



(x a ,f), y l = i(el), inE, 
(x a ,£,i), &=i,(ei), inE*, 



(x a ,y\x b ,yi) 

(X ,£i,X ,£,j) 

{x a ,y\p b ,n J ) 



in TE, 
in TE*, 
in T*E, 
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It is well known (cf. 
vector bundles: 



1*E 



(x a ,^,p b ,^) mT*E*. 
]) that the cotangent bundles T*E and T*E* are examples of double 

T*7T „ „ T* 





Note that the concept of a double vector bundle goes back to J. Pradincs 35, 55] , see also [31 • In 
particular, all arrows correspond to vector bundle structures and all pairs of vertical and horizontal 
arrows are vector bundle morphisms. The double vector bundles have been recently characterized |15) 
in a simple way as two vector bundle structures whose Euler vector fields commute. The above double 
vector bundles are canonically isomorphic with the isomorphism 

(2.1) K T :T*E — >T*E* 

being simultaneously an anti-symplectomorphism (cf. [51 1221 IT?] ). In local coordinates, 1Z T is given by 

K T (x a ,y l ,Pb,Tr j ) = (x a ,ir u -p b ,y J ). 

This means that we can identify coordinates Ttj with £j , coordinates ip J with y 3 , and use the coordinates 
(x a ,y l ,pb,Cj) hi T*E and the coordinates (x a ,£i,p b ,y J ) in T*E*, in full agreement with (|2.ip . 

For the standard concept and theory of Lie algebroids we refer to the survey article [26j (see also 
[Til [27]). It is well known that Lie algebroid structures on a vector bundle E correspond to linear 
Poisson tensors on E* . A 2-contravariant tensor II on E* is called linear if the corresponding mapping 
II: T*E* — » TE* induced by contraction, II(^) = i v H, is a morphism of double vector bundles. 
One can equivalently say that the corresponding bracket of functions is closed on (hber-wise) linear 
functions. The commutative diagram 




describes a one-to-one correspondence between linear 2-contravariant tensors II on E* and morphisms 
e (covering the identity on E*) of the following double vector bundles (cf. [121 E]) : 



(2.2) 



J*E- 



TE* 




TM 



In local coordinates, every such e is of the form 

(2.3) s(x a ,y\p b ,^ = (x a ,^p b k (x)y k ,c^(x)y^ k +a^x)p a ) 

(summation convention is used) and it corresponds to the linear tensor 

IT - 4(izi)a% ® % + P \{x)d ti ® 3 xb - a?{x)d x a ® % . 
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The morphisms (|2.2p of double vector bundles covering the identity on E* has been called an algebroid 
in [17] , while a Lie algebroid has turned out to be an algebroids for which the tensor II e is a Poisson 
tensor. We can consider the adjoint tensor 11+ , i.e. the 2-contravariant tensor obtained from II £ by 
transposition: 

11+ = c£.(a;)£ fe % ® 8 (] + p\{x)d xi <g> % - ^{x)d ij <g> 

and the opposite tensor — II E . It is clear that 11+ and — II £ are linear. They correspond therefore to 
new algebroid structures: the adjoint algebroid structure s + and the opposite algebroid structure e. An 
algebroid we call a quasi-Lie algebroid if e + = e. 

The relation to the canonical definition of Lie algebroid is given by the following theorem (cf. 

Theorem 1. An algebroid structure (E,e) can be equivalently defined as a bilinear bracket [-,-] E on 
the space Sec(E) of sections of t: E — > M, together with vector bundle morphisms p, a : E — > TM (left 
anchor and right anchor), such that 

[fX,gY] e = f ■ p(X)(g)Y - g ■ <r(Y)(f)X + fg[X,Y} £ 

for f,g G C°°{M), X, Y G Sec(i?). The bracket and anchors are related to the bracket {tp,ip}Yi e = 
(IL, dip ® dip) in the algebra of functions on E* which is associated with the 2-contravariant tensor H e 
by the formulae 

i{[X,Y] e ) ={i(X),i(Y)}n c , 
n*(p(X)(f)) ={o(X),n*f}n e , 
n*(a(X)(f)) ={n*f,L(X)} Us . 

The algebroid (E, e) is a quasi-Lie algebroid if and only if the tensor IL e is skew- symmetric, and it is 
a Lie algebroid if and only if the tensor Tl £ is a Poisson tensor. 

Since the dual bundles of tte ■ T*E — > E and Tn : TE* — > TM are, respectively, te : TE — > E 
and Tr : TE — > TM, the dual to £ is a relation k = n e : TE — \>TE. It is a uniquely defined smooth 
submanifold k in TE x TE consisting of pairs (v, v') such that p(te(v')) = Tt(v) and 

(v,e{v*)) jT = {v',v*) TE 

for any v* G T* , ,->E, where (•, -)j T is the canonical pairing between TE and TE* , and (•, -) TE is the 
canonical pairing between TE and T*E. We will write k : v — [> v' instead of (v, v 1 ) G k. This relation 
can be put into the following diagram of " double vector bundle relations" 

(2.4) TE< TE 




The relation 



TE< 




TE 



TM 



is a vector bundle morphism of the second kind, i.e. it is represented by linear maps of the fiber TE 
over v G TM into the fibers T e E for all e G E such that p(e) = v. This is also the simplest example 
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of a morphism of Lie groupoids in the sense introduced and exploited by S. Zakrzewski [46 . To 
such relations we will refer therefore as to Zakrzewski morphisms. The expression of the Zakrzewski 
morphism (|2.4p . dual to e, in local coordinates reads 



(2.5) k : (x a , Y\ p b k (x)y k , — > (x a ,y\ <r b k (x)Y k , Y' + c> kl {x)y k Y l ) . 

It is easy to see that the relation k" 1 coincides with K s +. Thus n = k^ 1 for quasi-Lie algebroids. 

A canonical example of a mapping e in the case of E — TM is given by e = Em = cx~m ~ ^ ne inverse 
to the Tulczyjew isomorphism olm '■ TT*M — » T*TM |38j . The dual Zakrzewski morphism is in this 
case the well-known 'canonical flip' km '■ TTM — > TTM. Since cum is an isomorphism, km is a true 
map, in fact - an isomorphism of the corresponding two vector bundle structures as well. 

A C 1 -curve 7 : M — > E (or a C 1 -path 7 : [toj^i] ~~ * E) in an algebroid E we call admissible, if the 
tangent prolongation t(7) of its projection 7 = 707 coincides with its anchor: 

(2.6) t(7) = P (j(t)). 

A curve (path) in the canonical Lie algebroid TM is admissible if and only if it is a tangent prolongation 
of its projection on M. If we denote T hol E the subset of TE consisting of holonomic vectors, 

(2.7) T hol E = {veTE: Tt(v) = p(r E (v))} , 

then admissible curves (paths) in the algebroid E can be characterized as those curves (paths) whose 
tangent prolongations lay in T hol E. The set of holonomic vectors T hol E can be equivalently char- 
acterized as the subset in TE which is mapped via Tp : TE — > TTM to classical holonomic vectors 
T 2 M = {u 6 TTM : km{u) = u}, that justifies the name. In other words, 

T hoi E= (T^-^M). 

Note also that T hol E is canonically an affine bundle over E modelled on the vertical bundle VE C TE. 
In local coordinates, T hol E as submanifold in TE is characterized by the equations x a — p°;(x)y l , so 
(x a , y l , y 3 ) can serve as local coordinates in T hol E. It is easy to see that, for quasi-Lie algebroids, 
K(T hol E) = T hol E. 

Let now 7 : [to,ti] — > E be a path and ( : [to>*i] ~ ¥ V.E C TE be a vertical vector field along 7, 
t b(C(0) = ^ is well known that \/E ~ E 1 ©m B, so vertical vectors at e 6 can be canonically 
identified with vectors of the fibre E T i e \. Thus, the vertical vector field £ can be identified with a path 
C,e in covering 7. We can consider now the tangent prolongation X.(C,e) to get a vector field along 
Ce- The operation ( 1— > t(Cs) associates with any path £ in Vi? a path t(Ce)in T.E. For t> £ TE, in 
turn, the family k(v) defines a vector field over p~ 1 (Tt(v)). More precisely, for every e € p~ 1 (Tr(v)) 
there is a unique vector k(v) c £ T e i? such that k(v) £ € k(v). We get the following. 

Theorem 2. 1/7 : [to, ti] E is an admissible path in E, then every vertical vector field £ : [to, ii] ~^ 
Vi? aZoni? 7 defines canonically a vector field 5cj ■ [tojti] - * T_B a/ong 7 fry 

(2.8) «J c7 (t) = «(t(CB)(*)) 7( t) ■ 

In /oca? coordinates, with 7 (t) = (x a (t) , y l (t)) and ((t) = (x a (t), y l (t), 0, f l {t)), 



(2-9) S a (t) = /'(tjo^t))^^*)) + ^(t) + c k 3 (x{t))y\t)P{t) ) tffc ( 7 (t)) . 

In other words, in local coordinates in TE, 



(2.10) <5 c7 (t) = (V(t),y l (t), /''(t)^*)), ^{t)+c k J {x{t))y\t)P{t)^ 
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The vertical vector fields £ along 7 we will call vertical variations or vertical virtual displacements of 
7 and the vector fields <5^7 along 7 - admissible variations or admissible virtual displacements. Note 
that the space V(7) vertical variations of 7 is canonically an (infinite-dimensional) vector space. 

Remark. In [8j [33], analogs of the admissible variations ^7 have been obtained (in Lie algebroid 
context, of course) from tangent lifts of time-dependent sections of E. The tangent lifts of sections have 
natural generalizations for general algebroids (THE]. We prefer, however, to define the admissible 
variation S^j directly by means of the vertical variation £ and the relation k, as being more fundamental 
and conceptually closer to the standard concepts of variations. 

3 Lagrangian formalism for general algebroids 

The double vector bundle morphism (|2.2p can serve as geometric background for generalized Lagrangian 
formalisms. 

The Lagrangian L : E — > R defines two smooth maps: the Legendre mapping: Xl : E — > E* , 
Xl = te* o e o dL, which is covered by the Tulczyjew differential Al : E — > IE* , A^ = eo dL: 

(3.1) T*E 

dL 

E- 

The lagrangian function L defines therefore the phase dynamics V = Al(E) C TE* which can be 
understood as an implicit differential equation on E*, solutions of which are 'phase trajectories' of 
the system f3 : M — > E* and satisfy t(/5) (t) € L. An analog of the Euler-Lagrange equation for curves 
7 : R -> E is then 

{E L ): t(A L 07) = A L 07. 

The equation (El) simply means that Al o 7 is an admissible curve in TE*, thus it is the tangent 
prolongation of Al ° 7- In local coordinates, T has the parametrization by (x a , y k ) via Al in the form 
(cf. ((231) ) 

(fit Fit Fit \ 

X -,—( X ,y),pl( X )y\ C ^ X )yi — ( X ,y)+^(x) — (x >y )J 

and the equation (El), for 7(t) = (x a (t) , y l (t)) , reads 

,„ oX , . da;" a . . fc d / dL \ k dL <9£ 

(3-3) W= ■^- = /4(*)y fc , d ^^J=4W^ + (T ^ ) ^ 

in the full agreement with [211 [301 @S], if only one takes into account that, for Lie algebroids, 
a c - = pj. As one can see from (|3 . 3f) . the solutions are automatically admissible curves in E, i.e. 
P(l(t)) — ^(t j)(t)- As a curve in the canonical Lie algebroid TM is admissible if and only if it is 
a tangent prolongation of its projection on M, first-order differential equations for admissible curves 
(paths) in TM may be viewed as certain second-order differential equations for curves (paths) in M. 
This explains why, classically, the Euler-Lagrange equations are regarded as second-order equations. 

Remark. The Tulczyjew differential : TM — > TT*M with a given Lagrangian function L on the 
canonical Lie algebroid E = TM is sometimes called the time evolution operator K (see [T)), as the 
first ideas of this operator go back to a work by S. Kamimura [5U]. This operator has been studied 
by several authors in many variational contexts, however, without recognition of its direct relation to 
a (Lie) algebroid structure. We named this map after W. M. Tulczyjew, since our understanding is 
based on his ideas [39]. 

The time-dependent version of the above picture is the following. Consider the direct product 
E = E x TR of the algebroid E with the canonical (Lie) algebroid TR equipped with canonical 



TE* 



E* 
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coordinates (t,i). The corresponding algebroid morphism is clearly the product of e and the inverse 
of the Tulczyjew isomorphism cxr: 

(3.4) e= (e, a^) : T*E = J*E x TTR -> JE* x TT*R = JE* . 

The affine hyperbundle Ar. = {(t, 1) E TR} of TR is a Lie affgebroid in the terminology of [11)1 [Hi [13] . 
Similarly, the afhne hyperbundle E\ = E x Ar in 25 is an affgebroid (so £7 = 2?i x R understood as the 
product in fibers is canonically a special affgebroid in the terminology of [13). The morphism (|3.4|) 
can be reduced then to 

(e, n Au ) : 1*E X = T*E x T*Ar -> T2T x ^ R C T(£* x R) . 

Identifying Ar with R in an obvious way, we obtain a morphism of double affine bundles [13j 

(3.5) e = (e, ttr) :T*(£xR) = T*£ x T*R -> TE* xTl = T(E* x R) , 

where ttr : T*R -> TR is defined by 7f R (t, s) = (t, 1) G TR. 

Here, we view E — E xR canonically as a vector bundle f :fi = £xl-iM xR over Mxl (the 
pull-back bundle of E with respect to the projection M x R — > M) and E* xRas its dual £7*. The 
time-dependent analog of the diagram (|3.ip defining the Tulczyjew differential, for the time-dependent 
Lagrangian L : E x R — > R reads 

(3.6) T*(£xl) ^T(FxM) . 

t v. ....- ^ 

Ah T (E*X«) 

£xl — *- FxR 

Although there is a canonical identification Ar — R, the use of ^4r explains the definition of holonomic 
vectors in this case: since T M (Exi E ) = J hol ExA R , we assume T M (£xI) = T hol ExR. This is due 
to the fact that the time-dependent picture is, in fact, an affgebroid picture (see [34l [37l [T3l [19l [43] ) . 
In other words, A L : E xR-> TE* x TR ~ T(E* x R) and X L : E x R -> E* x R read 

(3.7) A L (e,t) = (A i i(e),(t,l)) ) A L (e,f) = (A L *(e),t) , 

where we put L*(e) = L(e,t) and we canonically identified TR with R x R. If now 7 is a curve in E, 
then the nonautonomous Euler-Lagrange equation reads 

(3.8) {ET): t(A L o 7 )=A L o 7 ,, 

where j(t) = (~f(t),t) is a natural extension of 7 to £ x R. The nonautonomous Euler-Lagrange 
equation in coordinates takes formally the same form (|3 . 3[) . but now with L depending on t. 

Example 1. There are many examples based on Lie algebroids, see for instance [71 [24l [TBI [29l [33] , 

(a) For instance, for the canonical Lie algebroid and the corresponding morphism - the inverse of the 
Tulczyjew isomorphism [38 

e = : T*TM -> TT*M , 
with y a — x a , we get the traditional Euler Lagrange equations 

dx a _ .„ d_ ( dL_\ _ dL_ 

(b) For a Lie algebroid which is just a Lie algebra with structure constants c^- with respect to a chosen 
basis, we get the Euler-Poincare equations 

dt\dyj) ljU dy k ' 
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(c) True Lie algebroid examples are usually obtained as reductions of standard Lagrangian systems 
on tangent bundles, like the reduction of the rigid body to a dynamics on so(3,R). Another example 
of this kind is a homogeneous sphere of radius r > 0, mass m, and inertia k 2 about any axis, moving 
on a horizontal table without friction (thus, is the table rotating or not makes no difference). In an 
obvious way, the system lives in fact on the Lie algebroid r : TR 2 x so(3,R) — > R 2 with product Lie 
algebroid structure. In standard coordinates the algebroid morphism 

e : T* (TR 2 x so(3, K)) -> T (T*M 2 x so(3, R)*) 

reads: 

(3-9) e (x,y,x,y,L),p x ,p v ,p ± ,py,p u ) = 

{x,y,p ± ,Py,p u ,x,y,p x ,p y ,LO' i p L02 -L) 2 pu> 3 ,UiPu 3 -wsp^tU^puj -WiPo, a ) . 



The pure kinetic Lagrangian 



L = im (± 2 + y 2 + fc 2 (w 2 + w| + wj 2 )) 



induces the "free" dynamics 

— (mx) = 0, —(my) = 0, ^-(mk 2 uj) = . 
at at at 

Later we will add nonholonomic constraints to this picture. 

The above examples are associated with Lie algebroids, but some "nonholonomic constraints" on 
Lie algebroids may lead to Lagrangian systems on quasi-Lie algebroids. This is related to quasi-Poisson 
brackets associated with nonholonomic constraints [281 04] . 



Example 2. (Algebroid of linear constraints) Consider an algebroid structure on a vector bundle 
E equipped with a Riemannian metric (•, -)e and a vector subbundle C of E. Let P : E — » C be the 
orthogonal projection. We can choose a local basis of orthonormal sections (e,) = (e a ,eA) of E such 
that (e Q ) is a basis of local sections of C. According to the d'Alembert principle 5L(t('y)(t)) £ C°, 
where C° C E* is the annihilator of C, which in our case (cf. (|3.3[) ) takes the form 

for certain functions /i^, the constrained dynamics is locally written as 

(MO, ^ = 0, ^ = rtW , |(«) _«S, W _ . 

If we deal with a Lagrangian of " mechanical type" 

L = \{f) 2 -V(x), 
then = y A = and the equations (|3.10p reduce to 

that can be viewed as the Euler-Lagrange equations of the algebroid associated with the orthogonal 
projection of the tensor II e onto C* according to the orthogonal decomposition E* = C° © C* . Of 
course, even when E is a Lie algebroid, if C is not a Lie subalgebroid, the projected tensor is not a 
Poisson tensor and we deal with mechanics on a general algebroid, in fact a quasi-Lie algebroid in this 
case, since the projected Poisson tensor remains skew-symmetric. 
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4 Variational calculus 

For a general algebroid structure e on the vector bundle r : E — > M and a smooth Lagrangian function 
L : E — > M we will define a version of a variational calculus as follows. Our (infinite-dimensional) 
manifold A4 will be the space of all C 1 -paths 7 : [io,ii] — * E in E. Of course, like in the standard 
variational calculus, by curves through the path 7 £ M we mean C 1 -maps 

h : [t ,ti] xl3 (t, s) 1 ^ h(i, s) e 2? 

such that h(t,0) = j(t). Thus, the tangent space T^Ai - the space of all possible variations of 7 - is 
represented by ^(i,0), i.e. by continuous paths 6~/ : [to»*i] — * T2? covering 7 - vector fields along 7. 
The admissible paths form a subset M which is a submanifold in A4 in a natural sense, since a path 7 
is admissible if and only if t(7) C T hol E. As easily seen (see also [33 ), a vector field 67 : [to, h] — > 23 
along an admissible path 7 belongs to T 7 7V if and only if o t(£ 7 ) is tangent to T hol E, i.e. 

(4.1) K S (t(5 7 )(f)) e TT hol E c TTi? , 

where ke ■ TTE — > TT25 is the canonical flip. 

Note that we use here 'infinite-dimensional manifold' structures in a very intuitive sense. However, 
we could have put rigorously a Banach manifold structure on A4, Af, etc, similarly as it has been done 
in [33]. On the other hand, because the Implicit Function Theorem will be not used, a less formal 
language is completely satisfactory for our purposes, so we will skip technical complications associated 
with the Banach manifold setting. 

The Lagrangian L defines a differentiable function (action functional) Wl ■ M- — ► K by 

(4.2) W L (f)= I' L{l{t))dt. 

Completely classically, the differential of the action dW^^), paired with the tangent vector Sj, gives 

(4.3) <*7,dWi( 7 )>= f 1 (5 7 (t),dL( 7 (t)))dt. 

J to 

Now, we will make use of the algebroid structure on E and we will reduce the differential AWl to a 
distribution T> over the submanifold N in A4 consisting of admissible paths. For an admissible path 
7 : [to, ti] — > E, the space 2?( 7 ) C T^A4 of this distribution is exactly the space of admissible variations 
(virtual displacements) ^7 as they were defined in (|2.8[) . i.e. 

(4.4) 2?( 7 ) = {S a : C G V( 7 )} . 

In this sense, the space V( 7 ) of vertical variations, which is geometrically well-understood as the space 
of sections of the vertical bundle \/E along 7, is a model space for the space 2?( 7 ) of admissible 
(mechanical) variations which does not have so nice geometrical description in general. The reader 
can easily check that in the case of the canonical Lie algebroid E = TM the admissible variations we 
have just introduced coincide with variations of tangent prolongations of paths in M (with not fixed 
end-points yet), as they are understood in Classical Mechanics. The geometrical meaning of these 
variations is usually not understood being hidden behind the 'obvious' Lie algebroid structure on TM. 

Let us consider now the differential dW^ being restricted to T>. Our aim is to show its special 
realization, very similar to the one present in the standard variational calculus of Analytical Mechanics. 
Of special interest are variations d^j coming from the set V( 7 )o of paths £ that vanish at the end-points, 
((to) = 0, £(ti) = 0. They form a submanifold T>o of V and analogs of the standard Euler-Lagrange 
equations are obtained as equations for critical points of (dWL)\v i l - e - f° r such 7 G N that dWi( 7 ) 
vanish on T>o{l)- Note however, that in contrast with what has been done in [33], being interested in 
the infinitesimal picture only, we do not care about global homotopies inside the manifold of admissible 
paths. In fact, our distribution is not tangent to M in general, so even "infinitesimal homotopies" go 
outside M in the case of a general algebroid. This is due to the following observation. 
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Theorem 3. The distribution T> is tangent to the submanifold Af of admissible paths if and only if 
the right and the left anchor coincide, p = a . and they induce a homomorphism of brackets: 

(4-5) p([X,Y] s ) = [p(X),p(Y)] vf , 

where [-,-] v f is the bracket of vector fields. In particular, T> C T7V if(E,e) is a Lie algebroid. 

Proof. It is a matter of easy calculations to show that, according to (|2.10p . the vector field 8c,l along 
j(t) = (x(t), y(t)) satisfies (|4.1j) if and only if 

^(t) (a) - pk) ( X (t)) + p(t)y\t) {^Lfi - M ff « _ c*.A {x{t) ) = . 

Since the above should be satisfied for any admissible 7 and for any given x(t) = x(to), we can take 
/■'(to), ^£~(*o) an d y(to) arbitrary. Hence we get p = a and 

dab J n u dp\ rknh _d P ) d P \ kb _ 

The latter can be rewritten in the form p([ei,ej] E ) = [pi^i), p(ej)] v f, whence (|4.5[) . □ 

Remark. One develops often a variational calculus introducing homotopies as " paths in path spaces" 
satisfying certain boundary conditions - this is exactly how the variational calculus on Lie algebroids 
has been developed in [33j . However, this approach is much more restrictive when passing to con- 
straints. Let us only mention the existence of singular paths in the theory of linear nonholonomic 
constraints. In this case no real variation of a singular path is possible, so the differential calculus does 
not make sense any longer. On the other hand, the standard Euler-Lagrange equations are obtained as 
critical points of the action - so in fact only "infinitesimal homotopies", i.e. admissible variations are 
used. In the Lie algebroid case, the admissible homotopies can be taken as integral curves of variations. 
M. Crainic and R. L. Fcrnandes have related homotopies of admissible paths to flows of the complete 
lifts of time-dependent sections of the Lie algebroid in their work [8] on integration of Lie algebroids. 
They did not mention the variational calculus, but this integration is actually finding a manifold G(E) 
(Lie groupoid) that allows to represent the variational calculi on the Lie algebroid E as reductions on 
standard variational calculus on TG(E). Let us also point out that, contrary to the approaches by 
M. Crainic - R. L. Fernandes and E. Martinez [HJ, [33] , we work in full generality and we do not assume 
at the beginning that admissible variations come from vertical variations vanishing at the end-points. 

Since calculating dWL{l) on T> according to (|4.3p . we can divide our path into a finite number of 
smaller parts if needed, we can assume for simplicity that the path 7 lies in a single coordinate chart 
(x a , y l ), so we can write j(t) = (x a (t) , y l (t)) . That our path is admissible means now that 

(4-6) = #&(*))»*(*). 

For an admissible variation ^7, with £(t) = f l (t)ei( , y(t)), we have then 

(J c7 (t),dL( 7 (t))> = [f k (t) ■ <( 7 (i)) • J^(7(*)) + (v l (t) ■ 4(7(t)) • f\t) + ^(t)) • ^(7(*))' 
= f k (t) (<( 7 (i)) • g.( 7 (t)) +y\t) .4( 2 (t)) - |£( 7 (t)) - ^(7(*))) + 

d^ 

Writing X L : E — » E* , \ L (x,y) = §^j{x, y)el, for the vertical derivative (Legendre map) associated 
with L, and the variation of the Lagrangian along 7: 

(4.7) M(t( 7 )(t)) = (^(7W) • j£(7(*)) + V\t) ■ 4(lW) • |£(7(t)) " ^^(7(*))) 4 , 
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we can write 

(4.8) (£ a (t),dL( 7 (i))) = ±({ E (t),\L(j(t))) + {( E (t), *i(t( 7 )(*))> • 

According to (|3.2p . it is clear that <5L(t(7)(i)) = if and only if the image of the path dL(7(t)) under 
e is admissible in TE* , i.e. if and only if 7 satisfies the Euler-Lagrange equations (|3 . 3[) . 

In a more explicit form the variation of the Lagrangian can be viewed as a map 5L : T hol E — > _E* 
which in coordinates reads, 

(4.9) SL(x,y,y) = 

(a?(x)§±(x,y) +y i <$ i (x)%k(x,y) -y^( x )^^(x,y) -y k g§^(x >y j) e{ . 

A geometrical description of the variation of the Lagrangian is as follows. If v G TE is a holonomic 
vector, v G T hol E, then, as easily seen, 

A L o t e , TA L : -> TE* 

are bundle maps over Al : E ^ E* and 5L(«) = Ax,(te(u)) — TAl(u) is a vertical vector in T TE r v \E*. 
As the vertical bundle \JE* C TE 1 * is canonically isomorphic to E* ®m E* by means of the vertical 
lift, we can identify 5L(v) with a vector 8L(v) = (6L(v))e* from the fibre of E* over t(te(u)) G M 
which, expressed in coordinates, is exactly (|4.9p . In other words, 

(4.10) 5L = ((A L o t e - T\ L )\ T hoi E ) E ^ = ((e o dL o r £ - T(t e , °eo dL))\ Jh oi E ) E , . 

Using the obvious pairing between VE and \JE* based on the fact that the fibers over e and e* , 
respectively, are canonically dual spaces if r(e) = 7r(e*), we can write (|4.8|) equivalently in the form 

(4.11) (* a (t),dL( 7 (i))> = A(C(t),dL( 7 (t))) + (C(i), *£(t(7)(*))> • 
Integrating (14. 8p (or (|4.1ip ) we get 

(S a ,dW L ( 1 ))= [ 1 (S a (t),dL( 7 (t)))dt = (C E (t),\ L (l(t)))\%+ [\c E (t),SL(t( 7 )(t)))dt 

J to J to 

(4.12) = C(i)(7(*))ll;+ Ac(t),&(t(7)(*)))<tt. 

•/to 

Now, if C G V( 7 ) , then (dW^( 7 ), <5 C 7> = *i(t( 7 )(t))) dt. If C G V( 7 ) , thenr(i)C(*) G V( 7 ) 

for any function r : [<o,ii] — > K, so dW^(7) vanishes on £>o(7) if and only if <5E(t(7)) = 0. 
We can summarize the above observations as follows. 

Theorem 4. By means of the variational calculus for a general algebroid one can define the velocities- 
momenta correspondence (Legendre map): Al : E — > E* and the variation of the Lagrangian SL : 
T hol E — » E* , such that the derivative of the action functional dWi,^) is represented by 

(dW L fr),6 a )=(C E (t),\Mt)))\%+ Au*),<5£(t(7)(i))>dt. 

J t 

Moreover, the formula \4. lUty defines the Tulczyjew differential A^ : E — > TE* associated with the 
Lagrangian L. An admissible path 7 (t) = (x(t),y(t)) in E satisfies (5L(t( 7 )) = if and only if dW l 
vanishes on Pq( 7 ) and if and only if ^ satisfies the Euler-Lagrange equations 



For a given admissible path 7 : [io,£i] — > E, the values p(fo) = Ai( 7 (io)) and p(<i) = Al( 7 (£i)) 
represent the initial and the /inaZ momenta, and ?j 7 (i) = <5£(t( 7 )(t)) - the external force that we have 
to apply to make the system moving along the path 7. A standard way to obtain the dynamics in 
Analytical Mechanics is to look for critical points of the action functional with respect to admissible 
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variations <5<;7 which vanish at the end-points. In this way we obtain the Euler-Lagrange equations 
(|3.3p for admissible curves in the form <$L(t(7)(f)) = 0. In a more general setting, one can view the 
force defining equation 

(4.13) *L(t( 7 )(t)) = 7fc(f) 

as a differential equation for 7 if the external force 7y 7 (i) is given. In many cases this force is defined 
in path-independent way as a time-dependent field of forces F:£xl-> E*, n(F(e,t)) — r(e), i.e. 
r) 7 (t) = F( 7 (t),t). 

There is no real difference when we admit time-dependent Lagrangians L : E x M — > R, so that the 
action reads 

(4.14) W L { 1 )= I 1 L(j(t),t)dt . 



to 



The formula (|4.12j) just takes the form 



(4.15) {dW L ( 7 ),S a ) = (f(t) ■ ^(7(*),*)) 



t Jt 



'<&(*)> «£(t(7)(*),*)>dt, 



where 



(4.16) *L(t( 7 )(t),t) = (ofW)) • ^(7(*),t) +»'(*) ' ^W*).*) - S^ (7(t),t) ) ei ' 

In local coordinates, 

(4.17) SL(x,y,t,y) = 

(o* (x)§±(x,y,t) + y% (x) fj, (x, y, t) - y*p? (x) (x, y, t) - y k (x, y)-^(x,y,t)) e{ . 

The geometrical picture is based on (13 ,6[) . Now, SL : T hol E x M E 1 * is defined as the map whose 
vertical lift is 

(4.18) SL = m v o SL = (A L o t e - TAL)| T ho (jB><R = (e o dL ot^ - T(rg. o e o dL))| T(l „ !£;xR 

and the standard Euler-Lagrange equation for the time-dependent Lagrangian with a presence of 
external forces takes the form 

SL(t^)(t),t)=r, y (t). 

Again, the equation SL(t('y)(t), t) = means that the image of the path dJj{-y{t), t) in T*(E x R) under 
e is an admissible path (tangent prolongation) in T(E* x R). 



5 Constraints 

In view of the just developed variational calculus on general algebroids we can introduce, in principle, 
two types of constraints: the configuration constraints which are put in the "bundle of velocities" E, 
i.e. constrains for paths in J\f, and the virtual displacement constraints put for variations, i.e. for fibers 
of the admissible distribution V. As the admissible variations are also related to paths in E, the latter 
constraints can be defined also via constraints in E that often leads to misunderstandings. In this way, 
a constrained submanifold in E (classically in E — TM) is sometimes referred to as a nonholonomic 
constraint. Note however that in general, speaking on a submanifold (in general - subset) of E as of a 
constraint does not make much sense before wc decide how the constrained submanifold produces true 
constraints in the variational calculus. To put some order in the subject, we will start with describing 
our understanding of constraints in the variational calculus for general algebroid that will motivate a 
description of constraints in the pure geometric setting. 
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Definition 1. A constraint in the variational calculus for a general algebroid is a subset C of the 
bundle T>. The corresponding (dynamical) configuration constraint is the subset C^f obtained from 
C by the projection to AT. The constrained variational calculus is the study of the differential of the 
action functional AWl restricted to C, or Co = C f] V . 

It seems that the true variational constraints in physics strongly depend on the actual system we 
work with. In theory however, the variational constraints are often derived from geometric constraints 
of different types in a more or less canonical way. A geometric constraint will be understood as a 
submanifold (more generally - a subset) S in E. Of course, as we have already mentioned (see also 
[40]), the submanifold (subset) S C E does not define a true variational constraint without additional 
specifications. There are at least two geometrically justified ways of deriving variational constraints 
out of S. According to the tradition (see the review article [6]), we will refer to them, respectively, as to 
vakonomic and nonholonomic constraints. In the vakonomic case we accept only admissible variations 
(virtual displacements) which are tangent to the constraint, while in the nonholonomic case we admit 
only vertical variations which are tangent to the constraint, i.e. which belong to V(S') = TSf] ME. 

Definition 2. 

• The vakonomic constraint associated with S C E is the variational constraint C vk (S) consisting 
of these admissible variations S^j which are tangent to 5, i.e. S^(t) £ TS. In particular, the 
admissible path 7 lies in S. 

• The nonholonomic constraint associated with S C E is the variational constraint C nh {S) consist- 
ing of admissible variations ^7 associated with vertical variations C, which are tangent to S. In 
other words, ((t) € TS (thus, £(i) € V(S)). In particular, the admissible path 7 lies in S. 

• A geometric constraint S C E we call holonomic, if the nonholonomic variational constraint 
associated with S implies the vakonomic constraint, i.e. C nh (S) C C vk (S). 

Note that the variational constraints associated with S can be very small or even empty, e.g. when 
there are no admissible paths in S. To avoid pathologies like that, certain additional integrability 
conditions can be introduced. A natural integrability condition we will use is p(S) C TSm, where 
Sm = t(S). It is assumed in the sequel that the geometric constraints are integrable. 

Remark. We should stress here the obvious fact that TS" is well defined in a general setting even if 
S is not a submanifold of E, since it makes sense to speak about smooth curves in E with values in 
S. Note that, just by definition, for the vakonomic constraint only the restriction of the Lagrangian 
function L to S plays the role in the variational problem. The latter is not the case for nonholonomic 
constraints, except for the holonomic case. We can say that holonomic constraints are those nonholo- 
nomic constraints for which only the restrictions of the Lagrange functions to S play the role in the 
corresponding variational problems. One can easily derive from the form of the lift (|2.10p that a linear 
(integrable) constraint 5, i.e. a vector subbundle S C E, is holonomic if and only if the algebroid 
bracket [■, -] £ is closed on sections of S, i.e. S a subalgebroid in E. 

1. Vakonomic constraints — variational approach. The variational problem depends now on 
studying the differential of the action functional on C vk (S). A naive but instructive approach is that 
the corresponding constrained Euler-Lagrange equations describe admissible paths 7 in S which are 
critical points of Wl relative to the generalized distribution CQ k (S) — C vk (S) P| T>q, i.e. such that 
dW L (j) vanishes on all ^7 S C^ fe (5)( 7 ): 

(5.1) (dW L (~f),5 a )= [\dL( 1 (t)),5 a (t))dt= [\{ E (t),6L(t( 1 )(t)))dt = 

Jt Jt 

for all vertical vector fields C along 7, with C(^o) = 0, C(^i) = 0) an( i such that ^7 is tangent to S. 
Of course, it is hard to decide how large is CQ k (S). This is related to the difficult questions of the 
existence of singular or abnormal paths, etc., which cannot be solved in the whole generality. Leaving 
these questions aside, we will reduce ourselves to natural and geometric sufficient conditions ensuring 
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that a given admissible path satisfies f)5 . 1 [) . Namely, let us observe that if <& is a function vanishing on 
S, then, as Sqj is tangent to S, 

(d$( 7 (t)),<y c7 (t)) = o. 

Thus, if 

(5.2) 5(L-^ k )(t( 1 )(t),t)^0, 7 (i)G5, 

for certain /Lt& : [to, *i] — * K and for certain functions $ vanishing on (e.g. defining S according to 
the Implicit Function Theorem), then, according to (|5.ip applied to L := L — /J.k$ k , 

(dW L (l),S a ) = j£ (dL(7(t)) > * c7 (t)>d* 
= j£ ((dL - Mi b(t)d$ fe ) ( 7 (*)) , * C 7(*)) d* - £ <&(*)> <S (L - M***) (t( 7 ) (*).*)) di = - 

so (|5.1|) is satisfied. Such 7 we will call a normal solution of the vakonomic variational problem associ- 
ated with S C E. In the above procedure we can take as well a time-dependent Lagrangian L satisfying 
(|5.2|) . The latter does not depend directly on how big is Cg fc (S') and it simply means that the image of 
the path d (L - /u fc $ fc ) [l(t),t) in T*(E x E) under e is an admissible path (tangent prolongation) in 
T(E* x M). Motivated by the tradition we will regard the equation (|5.2|l as vakonomically constrained 
Euler- Lagrange equation. There is a clear analog of the above procedure also for time-dependent 
constraints. The non-autonomous vakonomic Euler-Lagrange equation takes in coordinates the form 

(5.3) $ k (x,y) = 0, ^=P a k (x)y k , 

(5.4) &%(x,v,t)- 4j i*)y l fK*, v,t) - J£ v,t) = 

Mt)^(x,y) + n k (t) ^ t ^(x,y)-ci j (x)y i ^ r (x,y)-a](x)^(x,y)) 

and reduces to the classical one for the canonical Lie algebroid E = TM (see e.g. [5]). In the 
above form, however, the vakonomic Euler-Lagrange equation is not easily seen to depend only on the 
restriction of L to S. Below we present a geometric approach clarifying this question. On the other 
hand, this approach seems to be more transparent not only at this point. 

2. Vakonomic constraints — geometric approach. Of course, one can take (|5.3|) as the Euler- 
Lagrange equation for the vakonomic constraints without referring to the variational calculus. It has 
the advantage that we do not care about possibly complicate constrained admissible variations. There 
is a nice geometric interpretation of these equations. For simplicity we will reduce to the autonomous 
case, so that L does not depend on t. Let us recall first that, with any submanifold S in E and any 
function L : S — > R one can associate canonically a lagrangian submanifold Sl in T*E defined by 

Sl = {ot e £ T*i? :eeS and (a e , v e ) — dL(v e ) for every v e € T e S} . 

If S = E, then Sl = dL(E), i.e. Sl reduces to the image of dL. We have the following. 

Theorem 5. A curve 7 : K — > E satisfies the vakonomic Euler-Lagrange equations i5.2\) associated 
with an autonomous Lagrangian L : E — > K if and only if it is a projection to S of a curve 7* : K — * Sl 
whose image under £ : T*E — > TE* is admissible (is a tangent prolongation of a curve in E*). In 
particular, the vakonomic Euler-Lagrange equations depend on the restriction of the Lagrangian to the 
constraint only. 

Proof. If a curve 7 satisfies (|5.2p . then 7 is admissible and lies in 5. Moreover, the curve 7 *(t) = 
d£( 7 (*)) - ^ A (t)d^ fe ( 7 (t)) in T*^ lies in S L , projects on 7 and is mapped trough e to an admissible 
curve. 

Conversely, if a curve 7 *(£) has the above properties, then there are (ik(t) such that 7*(i) = 
dL(7(t))-^ fc (i)d$ /c (7(i)). Since admissibility of e(7*(i)) is is equivalent to S (dL - ^ fc (i)d$ fc ) ( 7 (i)) = 
(Theorem HJ, the theorem follows. □ 
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One can also think that the vakonomic Euler-Lagrange equations are not equations on curves in E but 
on curves in Sl- Then, we can just consider the projections of the solutions onto E. 
The corresponding diagram is the following 



(5.5) 



T*E d S L 
A 



.-v 



■TE* 



->E* 



where is the relation which is the inverse of the projection (tte)\s l '■ $l —> S and = eor^. Like 
in the non-constrained case, a curve 7 in S satisfies the vakonomic E-L equation if it is related via A £ 
to an admissible curve in TE*. 

From the above it should be clear that the phase space for the vakonomic constraint S is 
te* (e(TT~ 1 (S))) and the phase dynamic associated with the Lagrangian L is e(Sl)- There is an obvious 
version of the above picture in the non-autonomous case. 

Example 3. (Pontryagin's Maximum Principle) For an algebroid (E, e) over M consider the 
product algebroid Ejj = E x TU. Considering an optimal control problem in which the manifold U 
plays the role of set of control parameters and associated with: 

1. an integrable constraint S defined by means of a [/-dependent section / : M X U E of E by 
(e,v) € S e = /(r(e),Tjf(«)) and 

2. a Lagrangian function L : S — > M depending only on the base, L(e, v) — L(r(e), tm{v))- 

In local coordinates (x a , u a ,y l , u 13 ) in Ejj and the adapted coordinates (x, u, y, u 7 p x ,p u , £, n) in T*Eu, 
the product algebroid morphism ejj = (e,Sm) reads 

Eu(x,u,y,u,p x ,p u ,£,ir) = (x,u^,Tr,p b k (x)y k ,u,c k J (x)y l £, k + (x)p x a , p u ) . 

The Lagrangian submanifold C T*Ejj consists of points 



x,u, f(x,u),u, 



dL 

dx 



so the phase (implicit) dynamics is given by £;y (<!?£), which is the set of points 



x,u,£,0 ,p b k (x)f k (x,u) ,u,c k Jx)f{x,u)£, k +cf°;{x) 



and the vakonomic Euler-Lagrange equations read 



dL 

dx" 



df_ 

dx a 



(x, u) 



dL 

du 



df 
du 



(x,u) 



(5.6) 
(5.7) 

(5.8) 



^ = p b k (x)f^x,u), 



= 0. 



The equations (15. 6|) and (15. 7|) describe the phase dynamics on E* x U associated with the Hamiltonian 
H(x, u, £) = f l (x, u)^i — L(x, u) via the tensor II e - trivially extended from E* to E* x U . The equation 
(|5.8p in turn, is the equation for critical points of this Hamiltonian with respect to the control variable 
u. In the classical case E = TM, the equations (|5.7|) and (|5.8p read 



0. 



d Pa /0L (dL df»\ 

-dT = {dx-«~ Pb dx-« ' {X ' U) > (du^- pb dv?^ X ' u) 



We recognize the Pontryagin's Maximum Principle in its normal differential form. For Lie algebroids, 
this principle was first proposed in |31j . 
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3. Nonholonomic constraints — variational approach. A naive but instructive approach is to 
assume in this case that the constrained Euler-Lagrange equations describe admissible paths 7 in S 
which are critical points of Wl relative to the generalized distribution C^ h {S) 1 i.e. such that (|5.ip is 
satisfied for all 8gy G CQ h (S). Again, we will not discuss the problem how large is Cq H (S). Recall that 
5a G C nh {S) means that C{t) G V(5), where V(5) = TSf] ME is the vertical part of TS. If S is a 
submanifold and V(5 I ) has constant rank, then the annihilator (V(5))° C T*E\ S is a vector subbundlc 
(over S) in T*^. In this case the quotient bundle T*-Big/(V(5)) is canonically isomorphic to the 
bundle V*(5) - dual to V(5). The latter, viewed as a subbundle in pr\ : E ©Adr E — * E in an obvious 
way, is called the bundle of virtual displacements in [5j Section 8]. Of course, (V(S')) can be viewed in a 
similar way as a subbundle in E ®m E* — > In this interpretation, which we will generally use in the 
sequel, (V(5))° C E* {e) , e G S, is the annihilator of (V(5)) e C £ r(e) and (V*(S)) e = ^ (e) /(V(5))°. 
It is obvious that (|5.ip is satisfied for all ^7 G C$ h {S) if (and not only if, in general) 

(5.9) ai(t( 7 )(*)) e (V(S))° (t) . 

This equation we will view as the constrained nonholonomic Euler-Lagrange equation. Again, it is not 
exactly equivalent to the variational principle in general as it gives only a sufficient condition for a 
relative critical point of the action functional Wl- In local coordinates, if «3> fe are functions defining 
the constraint S via equations <fr k (x,y) = 0, then V(5)° is generated by ^fc- at points of S, so the 
constrained nonholonomic Euler-Lagrange equation reads 

(5.10) $ k (x,y)=0, *£.=fi[x)y i 

(5.11) TtW^v) ~ 4- Wgk*.v) " °$(z)l&(x,v) = M(t)$(x,v) . 

For E — TM this is exactly the Chetaev principle and for E being an arbitrary Lie algebroid the 
equations (|5.10|) , (|5.11|) coincide with the equations associated with nonlinear nonholonomic constraints 
considered in [51 [35] . 

Example 4. (rolling ball) Consider now the celebrated example of a ball rolling on a rotating table 
(cf. [2, 5j), more precisely, of a homogeneous sphere of radius r > 0, mass m, and inertia about any 
axis k 2 , moving without sliding on a horizontal table which rotates with constant angular velocity fl. 
Like in Example[TJ we can recognize that the system lives on the Lie algebroid r : TR 2 x so(3, R) — » R 2 
with product Lie algebroid structure and is ruled by the pure kinetic Lagrangian 

L = -m (x 2 + y 2 + k 2 (ujj +u 2 + ujI)) , 

this time however with the presence of nonholonomic constraints 

^(a;, y,a;,2/,u;) = x - ru 2 + Qy = , 
® 2 (x, y,x,y,u) = y + rui - fly = Q. 

According to (|5.10[) and (|5.11[) . we get the constrained nonholonomic Euler-Lagrange equation in the 
form 

x - ru 2 + Oy = Q, y + no x - fly = 0, ±(mx) = Hi, ^(my)=/i 2 , 
^(mk 2 ^) = r^ 2 , ^(mfc 2 w 2 ) = -rfii, ±(mk 2 cj 3 ) =0, 

that easily implies 

k 2 n . 

fc 2 ft _ (l 

^ r 2 + & 2 
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If S is a linear constraint, i.e. S is a vector subbundle in E, then (\/(S)) e can be identified with S T r e \ 
and (V(5))g with <SS e ) C E*,-,. In this case the constrained nonholonomic Euler-Lagrange equation 
(|5.9p takes the form 

(5.12) 5L(t( 7 )(t)) G S° im , 

which is exactly the d'Alembert's principle of virtual work. The d'Alembert's principle for Lie algebroids 
was first proposed in [7]. 

More generally, assume that S = A is an affine constraint, i.e. A is an affine subbundle in E. Then, 
\l(A) e can be canonically identified with the fiber v(A) T ( e ) C i? r ( e ) of a vector bundle v(A) which 
serves as a model vector bundle of A. Hence, (V(S'))g can be identified with v(^)r(e) C ^v( e ) and the 
constrained nonholonomic Euler-Lagrange equation reads 

(5.13) JL(t( 7 )(i)) e (v(A))° 



4. Affine nonholonomic constraints geometric approach. Let us assume that S = A is 
an affine subbundle in E (over Am) satisfying the integrability condition p(A) C TAm- In this case 
(V(A)) e is constant along fibers of A and it coincides with the v(A) T ( e ). Let v(A)° be the annihilator 
of v(A) which is a subbundle in E* (over Am)- Let i V (A) : v (^) ^ be the inclusion of v(A) in E, let 
*v(A) : ^LAfcf ~ * v (^)* be the dual of i v ^)> an d let 

: T(^ M ) - Tv(A)* 

its tangent prolongation. According to the integrability condition p(A) C TAm, the image e(T*E^) 
lies in T(E* Am ) and the corresponding diagram is the following 

(5.14) T*E {A 

dL 

Ed A 



■T(^Am) ^^Tv(A)' 



MA)' 



E 



A? 



T~v(A)* 



The space v(A)* is the phase space for the nonholonomic constraint A with A^ : A — > v(^4)*, A^ = 
*v(A) ° as * ne constrained Legendre map, and A^ : A — ► Tv(A)*, with A^ = Ti*,^ o e o dL, serves 
as the constrained Tulczyjew differential. The set A^(A) C Tv(A)* is the phase dynamics associated 
with the Lagrangian L. The nonholonomic Euler-Lagrange equation is described as follows. 

Theorem 6. A curve 7 : M — > A satisfies the nonholonomic Euler-Lagrange equation SL(t(j)(t)) G 
v(A)°/ t % i/ and on??/ i/ £/ie curve A^(7(t)) m Tv(A)* is admissible (is the tangent prolongation of a 

curve in m{A)* ). 

Proof. Consider local coordinates (x 1 ) = (x l ,x L ) on a open set U of M such that ^4m is determined 
by the constraint x L = 0. A local basis {e a } a= i,... jn _ r of sections of v(A) together with a section eo 
of A we can extend to local sections of E and complete them to a local basis of sections {eo, e a , e a } 
of the vector bundle E. Then, in coordinates (x 1 ,y ) = (x l , x L ,y° ,y a ,y a ) adapted to this bases, the 
local equations defining the constrained subbundle A as an affine subbundle of E over Am are x L = 0, 
y° = 1, y a — 0, so points of A have coordinates {x % , 0, 1, y a , 0). Note that integrability of the constraint 
A means that p L (x) = and p L a {x) = at points x G Am- 

Taking local coordinates (x l ,y a ) on v(A) we may write £ v (a) : v (^4) ^ £ as Wa)^, J/ a ) = 
(x i ,Q,0,y a ,0) and i* (A) (z%0,£o,£o,&*) = (at 4 , Co), so 
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For the adapted local coordinates (x i , x L , y°, y a , y a ,p i: p L , £ , £„, £ Q ) in T*E, the map e reduced to 
(T*i£)u takes values in T(I?u M ) (integrability) and reads 

e(x\ 0, 1, y\ 0, Pl , U) = (**,<>, U, pL(x\ (% e +p> (x i , 0), 0, {c° B {x\ 0)y e + c^ B (x i , 0))£ D +a I B (x i , 0) Pl ) . 
Therefore 

(5.15) Ti* {A) oe(x\Q,l,y a ,Q,pi,U) = 
(x\^p{{x\Q)y- + pl(x\0), (cf fc (^,0)y e + c°(x\0)) fc, + oft*'', 0)pj) 

and 

(a 4 , y Q ) = TC (A) o e (V, 0, 1, y a , 0, j^.0, 1, y a , 0), j^^', 0, 1, y«, 0) 
I*, 0, 1, y», 0), 0)y e + P V\ 0), {c°(x\0)y e + c° b (x\0)) jp^, 0, 1, y a , 0) 

+^(x)^ 7 (^,0,l,y a ,0) 

Therefore, locally, the nonholonomic Euler-Lagrange equations read: 

(5.16) x' = 0, y° = l, y a = 0, ^ = p> e {x\ 0)if + pg(s*, 0) 
(5-17) ^(x 4 , 0,1,^,0) = 

{eg (x l ,0)y e + eg (x* , 0)) (x* , 0, 1 , y a , 0) + o{ (x* , 0) (a* , 0, 1 , y a , 0) . 

On the other hand, (|5.17p means that 5Lb = for all b, i.e. <5£ £ v(A)°. □ 

In the case of a Lie algebroid and linear constraints A=v(A) covering the whole M, when we have 
in particular Am — M, <r\ = p\, the previous equations are precisely the nonholonomic equations 
obtained in [5] (see Equations 3.8). 



Again, there is an obvious version of the above picture for a time-dependent Lagrangian based on 
In the nonholonomic case, however, we cannot restrict the Lagrangian function to the constraint, except 
for the case which is, in fact, holonomic. 

5. Holonomic constraints. In the nonholonomic case we can restrict the Lagrangian L to the 
constraint S if the geometric constraint is holonomic. Note however, that this does not imply au- 
tomatically that the corresponding vakonomic and nonholonomic Euler-Lagrange equations are the 
same, since the equations are not precisely variational (they describe only sufficient conditions that 
the variational principle holds true) and they are obtained in different ways. On the other hand, in 
the linear case holonomicity means that the vector subbundle S is closed with respect to the algebroid 
bracket. Since the constraints are assumed to be integrable, for the canonical Lie algebroid E — TM 
this means, in turn, that S = TSm, so the constraints are holonomic in the classical sense. More 
generally, assume that A is an affine constraint, i.e. A is an affine subbundle in E. 

Theorem 7. An affine constraint A in a quasi-Lie algebroid E is holonomic if and only if the algebroid 
bracket of sections of A is a section ofv(A). 

Proof. Let us choose a basis of sections ei and the corresponding linear coordinates (x a ,y l ) in E such 
that A is locally defined by equations y l = 0, i > r + 1, and y r+1 = 1 and let 7(i) = (x(t),y(t)) be an 
admissible path in A. Then, y i (t) = for i > r + 1 and y r+1 = 1. Moreover, £ is a vertical variation of 
7, £(x(t), y(t)) = f l (t)d y i, if and only if f l = for i > r. In view of (|2.10p . 5c-f is tangent to A only if 

^(t)+4(x{t))y i (t)r(t)=0 
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for k > r. But for any k > r 

|(i)+^4{ I (% i (i)/i(i) = ^4 +1)j (i(t))f( f )+^ 4( X (tw(t)r(t). 

i.j j<r i.j<r 

As y l (i), /■' (i) for i, j < r are arbitrary, c*j- = for k > r and i < r + 1, j < r. Since (A = — c^, they 
vanish also for fc > r and for all < r + 1. This means that the bracket of local sections [ei,ej] £ , 
i,j < r + 1 belongs to the span of {ei, . . . , e r }, i.e. is a section of v(A). But sections of A are of the 
form e r +i + J2i< r ( Pi( x ) e ij so their brackets are sections of v(A). The converse is obvious. □ 

According to the terminology of [10l [Til HI] , one can say that affine holonomic constraints in a Lie 
algebroid are Lie affgebroids. A correct geometric description of time-dependent systems and other 
systems, based on the idea of Lie affgebroid was first proposed in [37l [34] and developed in [TO ] [Ti l [12] . 

6. Affine holonomic constraints geometric approach. 

If A is a holonomic affine constraint, then, using local coordinates as above, we can prove anal- 
ogously to the proof of Theorem [7] that c^ b (x l ,0) and CQ b (x l ,0) can be non-zero only for D indexing 
a section of v(A), symbolically D — d, and that <r^(a; 4 ,0) = 0. Now, using the local form (|5.15|) of 
Ti*, A \ o e, we conclude that Ti*^ o e vanishes on the annihilator of TA. Hence, Ti*^ o e defines a 
map e A : T* A — > T\/(A)* and the diagram (|5.14p reduces to the following: 

(5.18) TM 

A - 

This time, however, only the restriction of L to A does matter. The phase space is v(A)* , the phase 
dynamics is implicitly defined as A.f;(A) — e A o dL(A) C Tv(A)*, and the Euler-lagrange equation for 
a curve 7 in A reads 

A£ o 7 = t(Af 07). 

In local coordinates 

(5.19) i 4 = 0, y = l, y Q = 0, ^ = /^(xSOjy + pg(x<,0) 

(5-20) i(w^,V a )) = 

(c d eb (x\0)y* + c d ob (x\ 0)) ^{x\ y a ) + o\{x\ Q)^( x \y a ) . 

The above equations (canonically reduced to A) are exactly the Euler-Lagrange equations for a (Lie) 
affgebroid obtained in [T31 [TS] . One can say that Geometrical Mechanics on a (Lie) affgebroid is just 
Geometrical Mechanics on (Lie) algebroid with a holonomic affine constraint. 
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